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s.^ . Abstract 

' We obtain an exact solution of the Dirac equation in (2 + l)-dimensions in the presence of a 

V/^ . constant magnetic field normal to the plane together with a two-dimensional Dirac-oscillator po- 

, tential coupling. The solution space consists of a positive and negative energy solution, each of 

^'j^ \ which splits into two disconnected subspaces depending on the sign of an azimuthal quantum num- 

' ber, A; = 0, ±1, ±2, • • • and whether the cyclotron frequency is larger or smaller than the oscillator 

' frequency. The spinor wavefunction is written in terms of the associated Laguerre polynomials. 

For negative k, the relativistic energy spectrum is infinitely degenerate due to the fact that it is 
• 1—1 ■ 

, independent of k. We compare our results with already published work and point out the relevance 

' of these findings to a systematic formulation of the relativistic quantum Hall effect in a confining 

potential. 
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1 Introduction 



Due to recent technological advances in nano-fabrication there were a lot of interest in the study of 
low dimensional quantum systems such as quantum wells, quantum wires and quantum dots [IJ. In 
particular, there has been considerable amount of work in recent years on semiconductors confined 
structures, which finds applications in electronic and optoelectronic devices. The application of a 
magnetic field perpendicular to the hetero-structure plane quantizes the energy levels in the plane, 
drastically affecting the density of states giving rise to the famous quantum Hall effect (QHE) [2j. 
The latter remains among the most interesting phenomena observed in physics because of its link to 
different theories and subjects. 

The stationary state associated with the motion of electrons in a uniform magnetic field is a 
well-known textbook problem [3j. It results in a sequence of the quantized Landau energy levels and 
associated wavefunctions characterizing the dynamics in the two-dimensional (2D) plane normal to the 
applied magnetic field. This quantization has important consequences in condensed matter physics 
ranging from the classical de Hass-van Alphen effect in metals [3] to QHE in semiconductors [2]. 
The relativistic extension of these models turned out to be of great importance in the description 
of 2D quantum phenomena such as QHE in graphene O [6l [3 E] • In fact, several condensed matter 
phenomena point out to the existence of a (2 + l)-dimensional energy spectrum determined by the 
relativistic Dirac equation [9]. For very recent works, one may consult references \10 \ [T2 t 113 ^ I14j. 

Motivated by different investigations on the Dirac fermions in (2 + l)-dimensions, we give an exact 
solution of a problem that has been studied at various levels by researchers dealing with different 
physical phenomena. We have done so by considering a relativistic particle subjected to an external 
magnetic field as well as to a confining potential. By introducing a similarity transformation, we show 
that the system can be diagonalized in a simple way. Solving the eigenvalue equation, we end up 
accounting for the full space of the eigenfunctions that includes all cases related to different physical 
settings. More precisely, from the nature of the problem we get separate angular and radial solutions. 
The radial equation leads to a "kinetic balance" relation between the two-spinor components. In fact, 
depending on the range of values of three physical quantities, the full solution space splits into eight 
disconnected subspaces as summarized in Table 1 in Sec. 2.2. This allowed us to obtain various 
solutions and emphasis similarities to, and differences from, already published work elsewhere |15| . 

On the other hand, we give discussions of our results based on different physical settings. In 
fact, we show that for week and strong magnetic field there is a symmetry that allows us to go from 
positive to negative energy solutions (states and spectrum). This can be done by interchanging the 
confinement frequency w with the cyclotrons lVc and vise versa. This suggests defining an effective 
magnetic field that produces the effective quantized Landau levels. In both cases, there is a degeneracy 
of the Landau levels where each quantum number n is /c-times degenerate; in analogy with the non- 
relativistic case |16| . For the intermediate magnetic field case, it is underlined that the degeneracy 
is possible. Finally, we compare our findings with those in a very significant work by Villalba and 
Maggiolo [15]. The full rich space of solutions suggested enabled us to carry out a deeper analysis 
in relation to various physical quantities. For instance, we obtained, as expected in the absence of 
applied voltage, a null current density for both directions in the Cartesian representation. However, 
this is not the case in polar coordinate. In fact, we show that the radial current vanishes whereas 
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the angular component does not. It is dependent on various physical parameters in the problem. 
These results are summarized in Table 3 showing clearly the dependence of these values on the given 
subspace. This may offer an alternative approach for a systematic study and understanding of the 
anomalous QHE [71 [8]. Additionally, we discuss the non-relativistic limit of the problem. 

The paper is organized as follows. In section 2, we give the theoretical formulation of the problem 
where a similarity transformation is introduced to simplify the process for obtaining the solutions 
(spinor wavefunction and energy spectrum). We use the "kinetic balance" relation to obtain a sec- 
ond order differential equation for one of the two-spinor components. The second spinor component 
is obtained from this using the "kinetic balance" relation. The relativistic energy eigenvalues and 
corresponding spinor wavefunctions are obtained as elements in the eight subspaces of the full and 
complete Hilbert space. In section 3, we discuss the physical meaning of our results and their poten- 
tial application to QHE. To analyze the transport properties of the system we determine the current 
density in section 4 and the non-relativistic case in section 5. Finally, we conclude by discussing the 
main results and possible extension of our work. 



2 Formulation of the problem 

We start by formulating the problem in terms of our language ^17j. This is done by introducing a 
similarity transformation of the Dirac equation in polar coordinates. This will be convenient to handle 
the "kinetic balance" relation and therefore derive the full spectrum as a complete Hilbert space. 

2.1 Hamiltonian system 

The problem of a charged particle moving in a constant magnetic field B = B z is a 2D problem in 
the plane normal to the field [the Cartesian (x,?/)-plane or cylindrical (r, 0)-plane]. In the relativistic 
units, h = c = 1, the Dirac equation in (2 + l)-dimensions for a spinor of charge e and mass m in the 
electromagnetic potential Afj_ = (Aq , A) reads as follows 

[i7^(a^ + ieA^) - m] V = 0, A^ = 0,l,2 (1) 

where the summation convention over repeated indices is used, 'y^ = (7^,7) are three unimodular 
square matrices satisfying the anti-commutation relation: 

where Q is the metric of Minkowski space-time, which is equal to diag(-|- — — ). A minimal irreducible 
matrix representation that satisfies this relation is taken as 7*^ = a^, ^ = ia where {crj}^^^ are the 
2x2 hermitian Pauli spin matrices: 

"-("o)' 

Equation ([1]) could be rewritten as 

i^il) = (^-ia-V + e6f A + eAQ + mf3^ii (4) 
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where a and f) are the hermitian matrices: a = ia^a^fi = a^. We will see below that the symmetry 
of the problem is preserved even if we introduce an additional coupling to the 2D Dirac-oscillator 
potential. This coupling is introduced by the substitution V — > V + mujfp where to is the oscilla- 
tor frequency. For time independent potentials, the two-component spinor wavefunction ip{t, r, 9) is 
written as 

'il^{t,r,e) =e-''^ip{r,e) (5) 

and (jl]) becomes the energy eigenvalue wave equation (Ti. — e)ip = where e is the relativistic energy. 
The Dirac Hamiltonian Ti. is the 2x2 matrix operator: 

H = Ho + ia^a ■ fHr + i o^o ■ 6 He (6) 

where (f, 9) are the unit vectors in cylindrical coordinates and 

Ho = eAo + mas, 

Hr = —idr + eAr — imujrcjs, (7) 

He = — de + eAe. 
r 

For regular solutions of ([3]), square integrability (with respect to the measure = rdrd9) and the 
boundary conditions require that tp{r, 9) satisfies 

^i^{r,9)\r=o =0, ^{9 + 2t:) = ^{9). (8) 

r— ♦oo 

To simplify the construction of the solution, we look for a local 2x2 similarity transformation 
A(r, 9) that maps the cylindrical projection of the Pauli matrices (<? • r, a ■ 9) into their canonical 
Cartesian representation ((Ti,(T2), respectively That is 

Aa-fA-'^ = ai, Ka ■ 9 Ar^ = (T2. (9) 

A 2 X 2 matrix that satisfies this requirement is 

K{r,9) = \{r,9)e^"'-'^ (10) 

where A(r, ^) is a 1 x 1 real function and the exponential is a 2 x 2 unitary matrix. The Dirac 
Hamiltonian ([6]) gets mapped into 

H = AHhr^ = Ho- (T2Hr + aiHe (11) 

where different operators are given by 

Ho = Ho, 

Hr = —i — j-^ + ieAr — imujras, (12) 

He = (^e - y - ^fT3^ + eAg 
with Afc = dk\- Therefore, the 2x2 Dirac Hamiltonian becomes 

^ f m + eAo 9r-^ + jf+ ieAr - mtur - ^ {de - ^) + eAg\ 

^ \-9r + ieAr + mtur - " x) + -m + eAo J ' 

(13) 
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Thus, hermiticity of (|13p requires that 

A« = 0, ^-1 = (14) 

and fixes the exact form of the modulus of similarity transformation to be A(r, 9) = yjr. It is interesting 
to note that turns out to be the integration measure in 2D cyhndrical coordinates. We could have 
eliminated the A factor in the definition of A in (jlOp by proposing that the new spinor wavefunction 
X be replaced by ■^x(r, ^). In that case, the transformation matrix A becomes simply ea""^^, which 
is unitary. However, making the presentation as above gave us a good opportunity to show (in a 
different approach) why is it customarily to take the radial component of the wavefunction in 2D 
cylindrical coordinates to be proportional to Finally, we obtain the (2 + l)-dimensional Dirac 
equation {H — e) x = for a charged spinor in static electromagnetic potential as 



m + eAo — e i9r + ieAr — mojr — ^de + eA, 



-dr — ieAr — miLir — ^dg + eAg —m + cAq — e 




(15) 



where x± ^-re the components of the transformed wavefunction \x) = A 1-0). This equation will be 
solved by choosing an appropriate gauge to end up with the full Hilbert space. 

2.2 Eigenvalues and wavef unctions 

Now, we specialize to the case where a constant magnetic field of strength B is applied at right angles 
to the (r, 0)-plane, which is B = B z. Therefore, the electromagnetic potential has the time and space 
components: 

^0 = 0, A{r,9) = ^Bre. (16) 
Consequently, (jl5p becomes completely separable and we can write the spinor wavefunction as 

X±{r,e)=(t>±{r)T{9). (17) 

Thus, the angular component satisfies —i^ = where ^ is a real separation constant giving the 
function: 

T{e) = ^e*«^. (18) 



On the other hand, the boundary condition ip{9 + 2tt) = ip{9) requires that e*2'^^e-*'^37r ^ ^^ich, 
in turn, demand that e*^'^^ = — 1 giving the following quantum number: 

^ = ^K, K = ±1,±3,±5--- . (19) 

Consequently, the Dirac equation for the two-component radial spinor is reduced to 

= (20) 



where the physical constant G is given by G = m (uJc — uj) and uJc is the cyclotron frequency "^c = f^- 
Thus, the presence of the 2D Dirac-oscillator coupling did, in fact, maintain the symmetry of the 
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problem as stated below Q. Moreover, its introduction is equivalent to changing the magnetic field 
as eB — > eB — 2muj. As a result of the wave equation (j20p . the two spinor components satisfy the 
"kinetic balance" relation: 

<P±{r) (21) 



0=p(r) 



+ - + Gr 
ar r 



e lb m 

where e 7^ itm. Therefore, the solution of the problem with the top/bottom sign corresponds to 
the positive/negative energy solution. Using the "kinetic balance" relation (j2ip to eliminate one 
component in terms of the other in ()20p results in the following Schrodinger-like differential equation 
for each spinor component: 

[-§r. + + G'r' +[m'-e' + G ^ ± 1)] } Mr) = 0- (22) 

Again, we stress that this equation gives only one radial spinor component. One must choose either 
the top or bottom sign to obtain the component that corresponds to the positive or negative energy 
solutions, respectively. The second component is obtained by substituting this into the "kinetic 
balance" relation (j2ip . Nonetheless, we only need to find one solution (the positive- or negative- 
energy solution), because the other is obtained by a simple map. For example, the following map 
takes the positive energy solution into the negative energy solution: 

e — > —£, K — > —K, G — > —G, (j)± — > 0ip (23) 

which, in fact, is the CVT transformation. Here the charge conjugation C means that e — > —e and 
uj — > —ijj or the exchange of lo and lOc- It is easy to check that the above map ()23|) originates from 
the fact that the Dirac equation (j20l) is invariant under such transformation. Hence, we just need to 
solve for positive energies and use the above transformation to obtain the negative energy solutions. 
The total spinor wavefunction reads as follows 

Mr, Q) = ^ e^«^e- 5^^^ V(r) (24) 
where (/>(r) has two components, such as 

0=f**V (25) 



Equation ()22p looks like the non-relativistic oscillator problem with a certain parameter map of the 
frequency, angular momentum, and energy. For regular solutions of (p2]) . the bound states will be of 
the form 

<A± ~ z^e-'l''K{z) (26) 

where i-^rXz) is the associated Laguerre polynomial of order n = 0, 1,2, ••• and z = p^r^. The 
constants {/x, z/, p} are real and related to the physical parameters B, u and ^. Square integrability 
and the boundary conditions require that 2p > ^ and u > —1. 

Substituting the ansatz ([26|) into ([22]) and using the differential equation for the Laguerre poly- 
nomial shown in Appendix A, we obtain four equations. Three of them determine the parameters 
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{/X, zv, p} and one determines the energy spectrum. The first three are 

1 



2/x 



P 



--G, 

e > 
e <0. 



(27) 



For regular solutions of (j22p . the it sign in the expression for u corresponds to ib^ > 0. Now, the 
fourth equation gives the following (positive and negative) energy spectra: 



±m\ 1 + 



2\G\ 



2n+l± ^—^ + ^{s + s') 



(28) 



where s = sign (G) = igi and s' = sign (^). The sign of G depends on whether the oscillator frequency 
uj is larger or smaller than the cyclotron frequency Uc- To compare our work with frequently used 
notation in the literature, we can replace the quantum number ^ by + |, where A; = 0, ±1, ±2, • • • 
and — > — ^ imply that k — > —k — 1. In that case, one may write the energy spectrum as positive 
eigenvalues: 

= m^l + ^[2n + l + s + k{s + s')] (29) 

and negative ones: 

' ™ (30) 



'n,k 



m\ 1 + [2n + 1 + s' + k{s + s')] 



where s' = +1 for k = 0. It is interesting to note that for < the spectrum is infinitely degenerate 
because it is independent of However, for S^G > the degeneracy is finite and equal to n + k + 1. 
Substituting the wavefunction parameters given by (j27p into the ansatz (|26p gives for e > 0: 



as well as for e < 0: 



(f)_{r) = xl*^^2 1 e 



1 T-l 



L-^-^ (x2) , 



/fc > 
A; < 



/t > 
A; < 



(31) 



(32) 



where x = ryJ\G\ and are normalization constants that depend on the physical quantities lo and 
LOc- The lower component is obtained by substituting (jSip and (j32p into the "kinetic balance" relation 
(|2ip . Doing so while exploiting the differential and recursion properties of the Laguerre polynomials 
(see Appendix A) we obtain the following for e > 0: 



. xl^2l e 2'^ 



A 



n,k 



{s - l)(n - A;)L-'^-i (x^) - (s + l)(n + l)L;^r^ (x^) 



> 
k<0. 



(33) 
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On the other hand, repeating the same calculation for the upper component of the negative energy 
solution gives the function 



^ f A-+ [(1 + s){n + k + 1)L^ + (1 - s)(n + (x^)] , > 

\ <_^x[(l + s)L-'^-i(x2)-2L-'^(x2)], k<0 

which could have also been obtained by applying the CVT map to ([55]) . Thus, the structure 
of the whole Hilbert space solution consists of eight disconnected spaces that could be displayed in 
tabular form as shown in Table 1: 



Frequency 


LO > UJc 


U> > UJc 


U! < UJc 


UJ < UJc 


Energy 


e > 


e < 


e>0 


e < 


Azimuth 


k>0 k <0 


k>0 k <0 


k>0 k <0 


/c > /c < 



Tablel: Full space solution. 

Using the standard definition, we calculate all involved normalization constants in the above wave- 
functions. These are summarized in the Table 2, where as stated above s = sign(G) = \G\/G: 



Energy 


Azimuth 


Normalization 


e > 


A; > 






^{n + k)\ 






-1 


e > 


A; < 


Ki,k = y 




ii{n — k)\ 


l+(JT^{2(n + l) + (fc + l)(s-l)} 






e < 


/c > 




'„.{ 


An + k + 1)! [l + — ^ {2(n + 1) + k{s + 1)}! 




-1 


e < 


A; < 


^n,k = 




|7r(n A: 1)! 1 + , .^'^l , {2(n A:) + n(s + 1)} 

[ [ (-^n.fc-*^) 







Table2: Normalization in terms of different physical quantities. 



3 Discussions 



It is worthwhile investigating the basic features of some limits of our results and their interesting 
underlying properties. We consider three different cases corresponding to the relative strength of the 
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external magnetic field (cyclotron freqency) to the oscillator coupling (oscillator frequency). We also 
demonstrate the added value of our results as opposed of others in literature, in particular the classic 
work of Villalba and Maggiolo [15] . 



3.1 Energy spectrum properties 

To investigate the underlying symmetry of the system, one may study the properties of quantum 
numbers pairs {n,k). However, these may not provide simple hints on the ordering of the energy 
eigenvalues e^^, with the exception of two limiting cases: the weak and strong field. 

3.1.1 Weak field case 

Suppose that the cyclotron frequency is much smaller than the oscillator frequency. That is, uJc <^ 
G « —TULO or s = —1. Thus, one obtains the following positive 



<fcU.«c. ~ m-^/l + ^ [2n + k{s' - 1)] 



and negative energy spectrum 



-m\ll + — [2n+l + s' + k(s'- 1)1. 
*' m 



(35) 



(36) 



Consequently, for A; > (i.e. s' = +1) the energy spectrum is (semi-) infinitely degenerate since it 
becomes independent of k. Moreover, the two spectra are related as 



However, for < (s' = — 1) we obtain 



(37) 



e 



„,fcl<..,«c.-m^/l + — -(n-A;) 



(38) 



It is also interesting to note that for k >Q there exits a positive energy zero mode corresponding to 

P+ I 



m with the following spinor wavefunction: 





1 2 




I ^/muj re* 1 exp 


mujr 

2 


:) 



(39) 



where the normalization: 



^+ = <+ = ^2{^A:![l + 2^(A: + l)]}"\ 



(40) 



These results are in good agreement with those of Dirac fermions in the plane in the presence of a 
constant perpendicular magnetic field. 
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3.1.2 Strong field case 

Now, if the cyclotron frequency is much larger than the oscillator frequency then G ^ muJc and we 
obtain the positive relativistic energy spectrum: 

<fel-c»a; ^ m^^l+'^-^[2{n + l) + k{l+s')] (41) 

as well as the negative one: 



- -my 1 + ^ [2n + (/c + 1)(1 + s% (42) 
They are related to each other as 

^n,feUc>u) = £ji_l_l^fe_lUc>u;- (43) 

Thus, in this case the infinite degeneracy of the spectrum corresponds to negative values of the 
azimuthal quantum number (i.e., s' = — 1) where e+ = e^^^. Here, a negative energy zero mode exits 
that corresponding to £q |a)c>a; = ~m with the following spinor wavefunction: 

n \ 

(44) 



^o(r,^) = ^ e-t'^^^ (V^re*^) ' exp -\ 



1 2 






— -mcocr 






2 







where the normalization: 



A- = A-- = ^2 {7r(-fc - 1)![1 - 2^k]y\ (45) 

Comparing the week and strong magnetic field limits one can conclude that the dominant frequency 
that controls the physics of the problem is interchanged between the oscillator and the magnetic field 
as to iOc- More precisely, the fc- independent infinitely degenerate energy spectra are related to each 
other as follows: 

la;c«;w,fe<0 = ^itfl l'^c>£^,fe>0) \oJc<^u>,k>0 = lwc>a;,fc<0 (46) 

where the quantum number n corresponds to the Landau level index. The existence of a zero-mode 
energy is now very clear for positive (negative) energy with A; > (A; < 0), respectively. 

3.1.3 Fine tuned case 

If the oscillator frequency is tuned to resonate with the cyclotron frequency ( LOc) then 

G = mA, where A = Uc — ui such that |A| ^ m. In this case, the relativistic energy spectrum 
approaches the nonrelativistic energy limit 

E = —(e'^-m'^) (47) 
2m^ ^ ^ ^ 



giving the quantity: 



s + s' , s — s' 



2n + IH — ± — h k{s + s') 



(48) 
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In this case, the energy spectrum degeneracy occurs whence the quantum numbers associated with 
the two states ipi and V'2 satisfy the relation 

k2-h~ 2 • ^ ^ 

That is, when the ratio of the shift in the principal quantum number is matched with the shift in the 
azimuthal number either up or down depending on the relative strength of the two frequencies and 
sign of k. 

3.2 Comparisons with other studies 

We compare our results with those in very similar studies found elsewhere in the literature. Such 
as the classic work by Villalba and Maggiolo |15j and, in particular, the energy spectrum and spinor 
wavefunction. As for the latter, our is identified with theirs according to 



In what follows, we summarize our remarks regarding few points in [15]: 




(50) 



1. The imaginary i is missing from some of the off diagonal entries in the Dirac equation (11), 
however, it was later corrected in (28) and (29). 

2. The relative strength of the cyclotron frequency iVc to the oscillator frequency u; (i.e., whether 
uJc is greater than or less than uj) is ignored. 

3. In addition, the negative energy solutions was also ignored altogether. Thus, only one fourth 
of the regular solution space, which consists of eight subspaces and whose structure is shown in 
Table 1, was obtained in [15j. The authors obtained only the two subspaces corresponding to 
e > and cUc > uj. 

4. The alternative signs in (31), which correspond to the sign of the energy, was confused with the 
independent signs for the physical parameter fj, (k in our notation). 

5. The relative number of nodes of the top to the bottom spinor components for < as given by 
equations ([55]) and ([Mj) is incompatible with the "kinetic balance" relation ([21]). 



4 Density of current 

We examine the behavior of the present system by analyzing the electric current density. Indeed, from 
our findings we can show 

(a) =i((T3<T). (51) 

For this calculation, we use the spinor wavefunction obtained above. This gives a null value in Cartesian 
coordinates, which is Jx = Jy = 0. This, of course, is expected since there is no net charge drift. As 
a reassuring exercise, we calculate the same current in cylindrical coordinates 

Jr = J • f = i{a^a • f), Jg = J • = i{a^a ■ 9). (52) 
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In this calculation, we employ the similarity transformation pO|) . The calculation gives = 0; how- 
ever Jg does not vanish having the components given in Table 3. This is due to the fact that the 
physical problem in cylindrical coordinates is for a charged particle confined to a circular motion due 
to the constant magnetic field. 



Energy 


Azimuth 


Angular Current Component 


e > 


fc > 


^r = ^ [n+(^ + l)^]{l + ^-^[2(n + fc + l)+n( 

n,k ' ^ n.k ' ' 




e>0 


A; < 




r 


e <0 


A; > 




e <0 


A; < 




-i 



TableS: Density of current for 4 subspaces. 

One could make a different analysis in terms of the physical quantities corresponding to different signs 
{s = lb) and for all the eight different subspaces. All these analysis could be used to give an interesting 
description for the anomalous QHE. 



5 Non-relativistic limit 



It is interesting to study the non-relativistic limit of our work to reproduce results already known in 
literature. This can be achieved by taking the limit m — > oo in the above findings. Now in the units 
h = c = 1, the non-relativistic problem has already been worked (see, for example, |16j): 

^2 



2m 



•1 o , 2 2 

Z o 



(53) 



The wavefunctions are (where w = and s = -|-1) 

1 ,' 



^n,a(r,0) = (-1)" X 



V^k y {n + \a\y. 



exp 



2/g 



ia0 



(54) 



where n = 0, 1, 2, • • • is the principal quantum number, a = 0, ±1, ±2, • • • the angular moment quan- 
tum number and /q = magnetic length. The corresponding energy eigenvalues are given 

by 

/ I /n/l -I- 1 \ ( ,) . 

(55) 



En,a = il\n^ ^ — I + —a. 



where is the frequency 17 = -^w^ -f- 40;^ . 
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To compare with the non-relativistic Hmit of our work, we take the Umit m — > oo and use the 
well-known non-relativistic energy formula E = [e^ — m^c?) /2m giving 




A; > 

n-h-tfi, A;<0. 



i±i . (56) 



The lb sign for A; < is a remnant of the positive/negative energy spectrum of the relativistic theory 
that is exhibited as a zero energy mode in the infinitely degenerate part of the spectrum. 



6 Conclusion 

The present paper was devoted to give a complete solution to the confined Dirac fermion system in 
the presence of a perpendicular magnetic field. Indeed, using a similarity transformation we have 
formulated our problem in terms of the polar coordinate representation that allows us to handle 
easily the "kinetic balance" relation between the two spinor components. One spinor component was 
obtained by solving a second order differential equation while the other component was obtained using 
the kinetic "balance" relation. It resulted in a full solution space made of 8 subspaces, which suggests 
that it is necessary to include all components of this subspace in the computations of any physical 
quantity. A failure to do so will result in erroneous conclusions. 

Our results were employed to discuss few important limiting cases: the weak, strong and fine tuned 
magnetic field cases. In particular, we showed that there is a symmetry between the negative and 
positive energy solutions. In the weak magnetic field case the system was shown to behave like a two 
dimensional Dirac system in the presence of an effective magnetic field controlled by the oscillator 
frequency w. To support our analysis we compared our findings favorable with those available in the 
literature and underlined the reason behind some of our differences. 

On the other hand, we analyzed the transport properties of the present system in terms of the 
current density. As expected, we found a null current in Cartesian coordinates, however, in polar 
coordinates the angular component of the current was non-vanishing. Finally, we studied the non- 
relativistic limit where known results were recovered. 

The emergence of the quantum Hall effect in graphene [3 opened a good opportunity not only 
for experimentalists but also for theorists as well. Because of the relativistic nature of the fermions in 
grapheme and due to some additional constraints, the appropriate mathematical system seem to be the 
massless Dirac fermions. However, our present work suggests that present theoretical investigations 
in the literature did not include adequately contributions from all solution parameter space [13] and 
hence will lead to incomplete, and sometime erroneous, results. Extending our present analysis to the 
massless Dirac fermion system will be desirable to put the theoretical approach to grapheme systems 
on firm grounds. 

Finally, we think it will be appropriate to look for the irregular solutions of the present problem. 
The importance of this issue comes from the fact that it will help us to construct the two-point Green 
function which is very much need in the calculation of many physical quantities and will enable us 
build the corresponding conformal theory. 
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Appendix A: Properties of the associated Laguerre polynomials 



The following are useful formulae and relations satisfied by the generalized orthogonal Laguerre poly- 
nomials Ll^{x) that are relevant to the developments carried out in this work. They are found in most 
textbooks on orthogonal polynomials [19]. We list them here for ease of reference. 
The differential equation: 



■ (i2 d 

x—-^ + {i'+l-x) — + n 
ax^ ax 



Ll{x) = Q 



(57) 



where x > 0, > —1 and n = 0, 1,2,---. They could be expressed in terms of the confluent 
hypergeometric function as 



r(n + i)r(i/ + 1; 

The associated three-term recursion relation: 



■i-Fi(-n; V + l;x). 



xLl = {2n + u + 1)L5; - (n + z.)L^_i - (n + l)L^+i 
Other useful recurrence relations: 

xK = {n + u)L--^-{n + l)L''-\ 



The differential formula: 



The orthogonality relation: 



TV _ TV + l TU+l 



{x)Ln{x)Lf^,{x)dx = ^ 



r(n + i) 



where 



p'^{x) = x" e 



(58) 
(59) 

(60) 
(61) 

(62) 

(63) 
(64) 
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